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Abstract 

Starting from the quasiclassical equations for non-equilibrium Green's func- 
tions we derive a simple kinetic equation that governs ac Josephson effect in 
a superconducting quantum point contact at small bias voltages. In contrast 
to existing approaches the kinetic equation is valid for voltages with arbitrary 
time dependence. We use this equation to calculate frequency-dependent lin- 
ear conductance, and dc I — V characteristics with and without microwave 
radiation for resistively shunted quantum point contacts. A novel feature of 
the I — V characteristics is the excess current 2I c /ir appearing at small volt- 
ages. An important by-product of our derivation is the analytical proof that 
the microscopic expression for the current coincides at arbitrary voltages with 
the expression that follows from the Bogolyubov-de Gennes equations, if one 
uses appropriate amplitudes of Andreev reflection which contain information 

about microscopic structure of the superconductors. 
PACS numbers: 74.50. +r, 74.80.Fp, 73.20 Dx 
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Point contacts between normal metals have simple Ohmic I— V characteristics regardless 
of their electron transparency D. In contrast to this, I — V characteristics of the super- 
conducting point contacts may be highly nonlinear even in the simplest situation of short 
constriction between two ideal BCS superconductors, and exhibit a non-trivial dependence 
on D. The origin of this complexity is the oscillating Josephson current, which makes elec- 
tron motion in the contact essentially inelastic Jl],|2j . Recently, there has been a considerable 
progress in calculation of both dc and ac f||-|5| components of current in such contacts. How- 
ever, the results were limited to the situation when the contact is biased with a constant (in 
time) voltage supplied by an ideal source with vanishing impedance. 

It is of interest to generalize the theory of electron transport in superconducting point 
contacts to the case of finite impedance of the voltage source as well as to time- dependent 
voltages. This generalization is particularly important in view of the fact that most experi- 
mental realizations of superconducting quantum point contacts ||[7| are based on the super- 
conductor/semiconductor heterojunctions which typically have relatively large impedance. 
Below we develop such a generalization which is valid for small bias voltages, V <C A/e, 
where A is the superconductor energy gap in the electrodes. 

We consider a ballistic quantum point contact with characteristic dimensions much 
smaller than both the elastic scattering length and coherence length of the superconducting 
electrodes. DC supercurrent in such a contact is known to be carried by the two discrete 
energy states with energies e± = ±Acosy?/2 inside the energy gap |8|||, where ip is the 
Josephson phase difference across the contact. These states are spatially localized in the 
contact region because of the Andreev reflection. At low voltages V <C A/e, dynamics of 
these states is slow on the frequency scale given by the energy gap, tfi = 2eV/h -C A/ft, and 
one could expect the ac Josephson effect in this regime to be described in terms of the same 
two quasistationary states. However, in contrast to the stationary regime (V = 0) when the 
occupation of these states is given simply by the equilibrium Fermi-Dirac probabilities, in 
the non-stationary situation the occupation of these states is quite non-trivial. 

We first discuss our final result, the kinetic equation which governs the evolution of 
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occupation probabilities p± of the levels e±. (The systematic development leading to this 
equation is presented in the last part of the paper.) Because of the normalization condition 
Yj±P± — 1) h is convenient to write the kinetic equation in terms of the difference of the 
two probabilities, p((p(t)) = p~ — p+- Kinetic equation for p(<p) is: 

P(^))=7(£)[n(£)-^(i))] 1 (1) 

where n(e) = tanh(e/2T) is the equilibrium value of p(e(ip)); 7(e) is the rate of quasi- 

particle exchange between the bulk electrodes and discreet levels in the constriction, and 

e = e(ip) = Acos<^/2. The rate 7 is roughly proportional to the subgap density of states 

in the superconducting electrodes; it vanishes in the ideal BCS case; if the gap is slightly 

smeared by finite electron-phonon interaction, 7 is given by the following expression |10| , [lT| : 

r 9(e /2 -A 2 ) (e - e') 3 cosh(e/2T) 
7(e) -ajde ^_ _ ^ _ cogh(e y 2T) • ( ) 

Here a is a constant determined by the parameters of electron-phonon interaction. To 
the kinetic equation ([]]) we should add a "boundary condition" which states that the level 
occupation reaches equilibrium as soon as the levels hit the gap edges, e = ±A (Fig. 1), 
that is: 

p(cp) = (-l) m n(A), for ip = 2vrm , m = 0,±1,... (3) 

We can take into account small reflection coefficient TZ of the point contact, 1, by 
including in the kinetic equation the Zener transitions between the two levels which occur 
at the point (p = 7rmod(27r) with the probability A ||. For vanishing external resistance the 
transition probability is A = exp{— 27iTZA/h | ip |}, where ip is taken at the transition point. 
Account of the Zener transitions is achieved by imposing one more boundary condition on 
p{if): 

p(ip + Osign(0)) = (2A — l)p(if — Osign(0)) , for <p = 7rmod(27r) . (4) 

The function p(ip(t)) given by the eqs. ([!]) - ([J) determines the current J(t) in the point 
contact: 



/(*) = ^sin^W)), (5) 

where Rn = nfi/Ne 2 is the normal-state resistance of the contact, and N is the number of 
transverse modes which are all assumed to be identical. 

Equations ([I]) - (||) allow us to describe dynamics of the point contact under arbitrary 
bias conditions. As a first example, we consider the linear response of the voltage-biased 
point contact to small oscillations of the Josephson phase difference around some stationary 
point (p , i.e. <p(t) = ip + <Puie~ tU3t , | cp u |<^C 1. Equations ([l]), (|5|) with this ip{t) give that the 
current oscillates around the stationary value fig] : 



hi^Po) = sin — tanh(-^-) , (6) 

^ ' eR N 2 V 2T ; ' K ' 

so that I = I s + l^e^ wt , and the frequency dependent linear conductance is: 

v( \ - Iuj - 2eiIlJJ 2ttA A sin 2 (y /2) jjep) i_ (po i/_£o_^ / 7 n 
^ J " K " ~ ^iv L 4Tcosh 2 (eo/2r)a; + i7(eo) 2 C ° S 2 tan ^2T jJ - 10 

where e = e(^o)- This equation generalizes the corresponding expression obtained by 
Zaitsev []J for large temperatures T ;§> A. (Note that eq. (0), as well as the kinetic equation 
(|I|) give only the leading terms in small relaxation rate 7.) 

We also can define the linear response to small dc voltage, V <C Kj/e. In this case the 
phase increases indefinitely, cp(t) = (fo + 2eVt/h, but deviation of the occupation probability 
p from equilibrium is still small. For such an evolution of ip eqs. ([!]) and @ give: 

at) = iMt)) + ^ 2 , 7(e)rcos ^ (e/2r) • W 



where e = e(<p(t)). This equation is the generalization of the recent result [13] to energy- 
dependent relaxation rate 7(e). We see that both types of linear response are sensitive 
functions of 7(e) and therefore they may be used to measure the subgap density of states in 
the superconductors. 

Non-linear response of the ballistic (TZ— 0) point contact to constant voltage V at zero 
temperature can also be obtained directly from eqs. (ffl), (El). At T = we obtain: 



mil, < Ifi < 7T , 

/(^(t)) = / c signysin^ (9) 

- 1 2exp(-j£ 7 (e)df)-l, vr<^<2vr, 



where to is the time at which ip(t) = it. Equation (||) shows that the effect of finite relaxation 
rate 7 is very similar to the effect of finite reflection probability 7Z (see eq. (11) in |J). 

As another application of the kinetic equation (|1]) we consider resistively shunted super- 
conducting quantum point contact biased by an external current I e (see inset in Fig. 3). For 
such a bias condition, the evolution equation for (p reads: 

^ = '•-'<*'»■ do) 

where current I(ip(t)) should be calculated from eqs. ([!]) - (H) self-consistently with <p(t), and 
R e is the shunting Ohmic resistance, which adiabatic approximation requires to be small: 
R e <C Rjy. We limit ourselves to low temperatures, T C A. 

For very small external resistances R e /R^ <C {(fry/A) ,71}, the rate of <p evolution is 
small and the current I(ip) is given by the stationary relation ([]). Applying this relation in 
eq. (|i0|), we conclude that the dc I — V characteristic of the point contact is given by the 
same relation as in the standard RSJ model (see, e.g., Sec. 4.2): 



(4 2 + (^) 2 ) 1/2 -^- (11) 

tie, tie. 



In the opposite limit of relatively large external resistances R e jR?q 3> {(Try/A) , 71}, the 
current-phase relation I(ip) is the same as with no relaxation [Q: 

J(^)=/ c sign(\/)|sin||, (12) 

and the I — V characteristic is given by the following relations: 

(/2 _ J2U/2 

v= [±ic) I = I e -V/R e . (13) 

Re 4 arctan J(I e + I c )/(I e - I c )} 



The main qualitative difference between expression ( |i~3"|) and the quasistationary expression 
(|TT|) is the excess current: / — > 2/irI c for V ^> I c R e in eq. flTS]), whereas in eq. ( |Tl"D the 
current obviously vanishes at large voltages. Note that eqs. (|TT|) and (|T3|) , and Figs. 2 and 



3 below give the I — V characteristics in the form (i.e., without the linear term V/R e ) 
that is directly applicable to typical bias conditions of point contacts which are not shunted 
intentionally. Under these conditions there is no shunting resistance at zero frequency, but 
there is finite impedance of the biasing leads in series with the contact at frequencies of the 
Josephson oscillations. It is known that this situation can be reduced to the RSJ model by 
simple subtraction of the dc current through the resistor (see, e.g., fl4jl , Sec. 12.4). 



Figure 2 shows how I — V characteristics evolve from eq. (|TT|) into eq. (|13|) with increasing 
external resistance. The curves were calculated numerically from eq. (||) and fllOD assuming 
no reflection in the point contact (7Z= 0), and also assuming that 7 is a phenomenological 
constant independent of energy. In the case when this transition is driven not by finite 
relaxation rate 7 but by finite reflection TZ the curves look qualitatively very similar. 

Figure 3 shows dc I — V characteristics of the point contact under microwave irradiation, 
I e (t) = Iq + A cos(fit), which exhibit the usual Shapiro steps at voltages Vfc >m = (k/m)hfl/2e. 



In the standard way (see Ref. |TJ]], Sec. 10.2) one can get that at large frequencies the I — V 
characteristic in the vicinity of the steps has a hyperbolic shape (similar to eq. QTTj)) with 
the step height I c \ a m Jk(2eR e Am/hQ) |. Here are the Bessel functions, and a m are the 
Fourier expansion coefficients of I{(f), 

8 (-l) m m 4 1 

or a r , 



7r 1 — 4m 2 ' 7r 4m 2 — 1 ' 

for the quasistationary current and for the current given by eq. (|i~2"D , respectively. From 
these expressions we can conclude that the height of the subharmonic steps (m/1) which 
are the hallmark of the presence of higher harmonics in I(ip) depends strongly on the value of 
external resistance. It increases for small external resistance due to the current discontinuity 
at (p = 7r mod27r. 

Now we briefly outline the major steps leading to our basic kinetic equation (|1|). We 
start with the quasiclassical equation for non-equilibrium Green's functions of the super- 
conductors (for general introduction to this technique see, e.g., ]T5|]). The Green's functions 
can be represented as + G, where G is a space-dependent non-equilibrium addition to 
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the equilibrium part G^ that is constant inside each electrode. For short constrictions, 
equations for the retarded and advanced parts of G read : 



iVF ~if~ = ^ Hr ' Ai Gra ^ ' Hr ' A = ( Sr ' A + ^^^RA > ( 14 ) 

where Sr jA = ((e ± ^7i) 2 — (A ± i^) 2 ) 1 ^ 2 ', 7i,2 and 7^ are, respectively, inelastic and elastic 
scattering rates, vp is the Fermi velocity, and coordinate z measures the distance from the 
point contact (z = 0) into the electrodes (z — > ±00). All functions in eq. (|T4]) are matrices in 
the electron-hole space; for instance, G^\(e,e') = [((e±ryi)cr 2 + (Ai^)^)/^^]^ 6- e'), 
with cr's here and below denoting Pauli matrices. 

The functions G should decay inside the electrodes (at z —>■ 00). If we perform "rotation" 
in the electron- hole space diagonalizing Gr A : 

G ( $ A (z, - U R , A (e)G^ = ±^ - , ^ = 1 i fl *f ' , (15) 

1 ~~ a _R,A 



where = ( e ±^7i _ ^r,a)/( A ±272), eq. fll4|) shows then explicitly that solutions decaying 
inside the electrodes should have the following matrix form: 

= U^u£%±U R , G^ 2) = UX^'tUa , (16) 

where a± = <j x ± ia y , and G^ 1 ' 2 ' denote the function in the first (z < 0) and the second 
(z > 0) electrode, respectively. 

The total Green's functions should be continuous at the point contact (z = 0). Imposing 
this condition and taking into account that there is a voltage drop V between the two 
electrodes of the point contact we can determine the functions u RA in eq. flip]). At small 
voltages, V —>■ 0, we get then for the total Green's functions Gr iA = Gr )A + Gr A at -2 = 0: 

G R A {e, t) = [ ^G RA (e + j, e - = ± a * C °^/ 2 ~ arcco ^/ A )J + ^ , (17) 

h, A \ > ) J 2n h, A \ 2 , 2 j 2?r sin ^/ 2 - arccos(e/A) ±i0] V ; 

where G# 5 a depend on time t via the time dependence of the Josephson phase difference <p, 
ip = 2eVt/h + (fo . In eq. ( ]T7| ) we neglected the relaxation rates 7^2 assuming that they 
are small. This is a legitimate approximation since, as usual, the effect of small energy 



relaxation on the occupation probabilities (i.e., on Gk) is much more important than the 
effect on the density of states. For G R: a given by eq. fllT]), the latter effect would be a small 
broadening of the Andreev-bound level. 

One can check directly from eq. (|T7| ) that this equation agrees with the stationary Green's 
functions calculated first by Kulik and Omel'yanchuk ]T2" |. In particular, in the subgap 



range | e |< A it corresponds precisely to one of the two Andreev-bound discrete energy 
levels: ReG^ oc S(e — A sign(sin y>/2) cosip/2). Since the evolution equation ([14]) and, 
consequently, eq. fll7|) refer to electrons moving in the positive ^-direction (p z > 0) this is 
the level that carries current in one direction. Evolution equation for electrons with p z < 
differs only by the sign in front of vp . In this case we get that G Rt A corresponds to the 
energy level at e = — A sign(sinyj/2) cos<^/2. 

To find the current I(t) in the point contact we need to calculate the Keldysh component 
Gk of the Green's function |T0|,p|: 



/(*) = a|^/ deSv{a z (G%* >0 \e,t) - G$* <0 \e,t)} . (18) 

Equation for G% is: 
f)G 

iVF—jr- = H R°k + H k G * - G R R k - GkH a , H K = H R n - nH A , (19) 

where H R a are defined in eq. fll4]), and n is the equilibrium quasiparticle distribution, 
n(e, e') = tanh(e/2T)5(e — e'). This equation shows that Gk can be written as Gk = 
Grti — tiGa + Gh, where Gh is the part that satisfies the homogeneous equation: 

BG 

iv F -^ = H R G H -G H H A , (20) 

Following the same steps that led to eq. (|IED we get that Gh should have the following 
matrix form: 

G^' 2) = U R 1 u^ 2 \l±a z )U A . (21) 

Imposing again the continuity condition at z = 0, we calculate u H ' and then find the 
current at arbitrary voltages: 
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I(t) = £ I ke ™^ , 
k 



eR 



N 



n+2k 



detanh(— )(1- | a R (e) | 2 ) E II I a ^ e + meV ) ? II a n( e + meV ) 

^ n=0 m=l m=n+l 

(22) 

Equation (^) has the same form as the corresponding expression that follows from 
calculations based on the Bogolyubov-de Gennes equations 0. The only difference is that 
the function or(c) which has the meaning of generalized Andreev reflection amplitude now 
contains full information about the microscopic properties of the superconducting electrodes, 
and is, in general, different from its "ideal" BCS value. In the particular case considered 



here it includes finite energy relaxation rates 7^2. For small 7's, part of the eq. fl22|) related 
to the dc current (k = 0) reduces to the so-called BTK expression for the current 0. To the 
best of our knowledge, this is the first explicit proof that the widely used BTK approach is 
equivalent to the microscopic theory of electron transport in short ballistic constrictions. 

Finally, to obtain the kinetic equation (ffj) we consider the Green's function Gk in the 
limit V — > 0, when it is given by the following expression: 

G K (e,t) = — I sin ^ I 5(e - A sign(sin<^/2) cosy9/2)iV(e, t) , 

where 

N [e , *)=»(«)+/ d^M-l eVVA T_> (23) 

Here ^7 = 2(71 — (e/A)7 2 ), so that 7 is given by the eq. @. 

Comparison of this expression for Gk with the subgap density of states that follows from 
eq. ( |iT|) shows directly that iV(e, t) has the meaning of quasiparticle distribution, so that 
(1 — N(e,t))/2 can be interpreted as an occupation probability of one of the two Andreev- 
bound levels inside the gap. Equation (|23|) with this interpretation immediately gives the 
kinetic eq. (fl|) and the boundary condition (||). Indeed, taking into account the definition 
of p in the kinetic equation we see that it is related to N as follows: p = A^sign(siny?/2). 
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This relation together with the eq. (p3|) give the boundary condition (^). Furthermore, 
differentiating eq. fl23|) with respect to energy and making use of the relation between energy 
and phase, e = sign(sin<£>/2)Acos<^/2, we get: 



If we express iV in this equation in terms of p we finally arrive at eq. ([!]). Although we 
assumed so far that the voltage V is constant in time, it is obvious that the evolution 
equations (|24|), (|l|) in the differential form are valid for arbitrary time dependence of the 
voltage, as long as the voltage itself and the rate of its variations are small. 

As a last remark we should mention that thermalization of the occupation probability p 
due to the boundary condition ([3D is instantaneous only on the long time scale set by the 
period of the Josephson oscillation. Crude estimate of the energy interval Se near the gap 
edge which determines p is (Ae 2 ^ 2 ) 1 / 3 , so that the corresponding time scale of thermalization 
is St ~ h/5e. In the relevant limit eV/A — > 0, 5t is much less than the period of the Josephson 
oscillations. 

In conclusion, we developed an adiabatic theory of the ac Josephson effect in short con- 
strictions between two superconductors. The theory is based on the simple kinetic equation 
for the non-equilibrium occupation probabilities of the two Andreev-bound states localized 
in the constriction. The kinetic equation is rigorously derived from the microscopic equa- 
tions for quasiclassical Green's functions of the constriction, and can be applied to situations 
with arbitrary time dependence of the bias voltage. 
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FIGURES 



Figure 1. Energies e± of the two Andreev-bound levels in a short constriction between 
two superconductors as functions of the Josephson phase difference <p. Solid dots represent 
"thermalization points" where occupation of the two level always reaches equilibrium. The 
diagram illustrates sign convention for the kinetic equation (J^) and the boundary condition 

Figure 2. DC I — V characteristics of the resistively shunted ballistic superconducting 
quantum point contact for various ratios of external Ohmic resistance R e and relaxation 
rate 7 at zero temperature; u c = 2eI c R e /h. 

Figure 3. Effect of microwave radiation with the amplitude A = 21 c and frequency 
Q = 2u c on the DC I — V characteristics of the resistively shunted superconducting point 
contact. The contact parameters are the same as in Fig. 2. From top to bottom, the curves 
correspond to "f/uj c = 0, 1, 00. 
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